Spin Correlations and Magnetic Order in Nonsuperconducting Nd2-xCea;Cu04±5 
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We report quantitative neutron scattering measurements of the evolution with doping of the 
Neel temperature, the antiferromagnetic correlations, and the ordered moment of as-grown, non- 
superconducting Nd2-i:Cei,Cu04±i (0 < a:: < 0.18). The instantaneous correlation length can be 
effectively described by our quantum Monte Carlo calculations for the randomly site-diluted nearest- 
neighbor spin- 1/2 square-lattice Heisenberg antiferromagnet. However, quantum fluctuations have a 
stronger effect on the ordered moment, which decreases more rapidly than for the quenched-disorder 
model. 

PACS numbers: 74.25.Ha, 75.40.Mg, 75.50.Ee 



Over the past eighteen years much has been learned 
about the evolution with doping of magnetic correlations 
in the high-temperature superconductors. For the proto- 
typical material La2-xSra;Cu04 (LSCO), neutron scat- 
tering has led to the observation of the loss of long-range 
antiferromagnetic order and the concomitant decrease 
of two-dimensional (2D) spin correlations .1], an incom- 
mensurate spin-density wave 2] , and coupled charge and 
magnetic "stripe" order upon codoping with Nd 3] . The 
richness of the observed behavior is attributable to the 
evolution of the cuprate superconductors from Mott- 
insulating parent materials which are found to be excel- 
lent realizations, at low energies and long wavelengths, 
of the spin-1/2 square-lattice Heisenberg antiferromag- 
net (SLHAF) 4]. Although the vast majority of research 
has focused on the hole-doped systems, the electron- 
doped materials, typified by Nd2-xCea;Cu04±5 (NCCO) 
0, have presented an important challenge to the field 
of high-Tc research and particularly to the paradigm es- 
tablished from experiments on LSCO. Unlike for LSCO, 
the magnetic phase of NCCO is very robust with respect 
to doping [a, l3,Ja and the magnetic correlations remain 
commensurate Isfis!. 

Since superconductivity in the high-Tc cuprates ap- 
pears in close proximity to antiferromagnetic phases, it 
is essential to understand the nature of nearby mag- 
netic ground states. As-grown, NCCO is not supercon- 
ducting, and the antiferromagnetic phase extends to the 
highest concentration for which samples can be produced 
(x « 0.18). Oxygen impurities are believed to occupy the 
nominally vacant apical sites |lfll | , and a relatively severe 
oxygen reduction procedure must be applied to induce 
superconductivity above x = 0.13 0. Although the re- 
duction procedure weakens the magnetism, magnetic or- 
der is still observed near optimal doping where — 24 
K 0, [l3| • One qualitative argument often employed to 
explain the robustness of the antiferromagnetic phase in 



the electron-doped cuprates is that, whereas hole carriers 
are doped into the O 2p band and frustrate the antiferro- 
magnetic arrangement of the neighboring Cu ions, doped 
electrons primarily reside on the Cu site, filling the 3d 
shell and removing the spin degree of freedom la 111 ■ 
New, quantitative insights have been gained recently into 
the role of static non- mag netic irnpurities introduced into 
the spin-1/2 SLHAF [i2ElElIli, and it would be in- 
teresting to test how well NCCO can be described by this 
model. 

In this Letter, we present detailed neutron scatter- 
ing results on the evolution with Ce doping of the 
magnetic properties of as-grown, non-superconducting 
Nd2-:rCe:rCu04±5 (0 < x < 0.18). Wc find that the in- 
stantaneous spin correlations in the paramagnetic phase 
are, in effect, well described by the randomly-diluted 
spin-1/2 nearest-neighbor (NN) SLHAF model for which 
we have carried out numerical simulations. On the other 
hand, the ordered moment falls off more rapidly, a clear 
indication that quantum fluctuations are stronger than 
for the model magnet and that such fluctuations mani- 
fest themselves differently for different observables. We 
also find that the primary difference between the Neel 
temperature and ordered moment of as-grown and re- 
duced NCCO is a shift in electron concentration that 
corresponds to |Aa;| « 0.03. 

Crystals of Nd2-2:Ce2:Cu04±5 were grown in 4 atm of 
O2 using the traveling-solvent floating-zone technique. 
The solubility limit is about x = 0.18, and our study 
spans the range < x < 0.18. In addition, several sam- 
ples were reduced for 20 h at 960°C in Ar, followed by an 
anneal for 20 h at 500°C in O2, for additional magnetic 
order parameter measurements. Magnetometry indicates 
that the reduced x = 0.14 and x = 0.18 samples are su- 
perconducting with Tc = 24 and 20 K (onset), respec- 
tively. Our single-grain crystals are cylindrical, typically 
with a diameter of 4.5 mm, and range in volume from 
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FIG. 1: Representative two-axis scans for x = 0.10 and 0.14 
along q2D = {Hl2 — 1/2, H/2 — 1/2) in the paramagnetic 
phase with collimations 40'-42.5'-sample-23.7' and incident 
neutron energy Ei = 14.7 meV. The instrumental resolution 
is indicated by the dashed lines. 



0.1 to 0.7 cm'^. Average Ce concentrations were verified 
using inductively coupled plasma spectroscopy and typ- 
ically found to equal the nominal values. The neutron 
scattering measurements were performed on the thermal 
triple-axis instruments BT2 and BT9 of the National In- 
stitute of Standards and Technology Center for Neutron 
Research. 

All as-grown samples exhibit long-range magnetic or- 
der at low temperature. However, significant two- 
dimensional (2D) magnetic correlations are already 
present well above the Neel temperature These may 
be probed through "two-axis" scans in which one inte- 
grates the dynamic structure factor S(c\2£,,uj) over the 
limits set by the thermal energy, —k^T, and the in- 
cident neutron energy, Ei. If the integration is car- 
ried out over the relevant fluctuations one obtains the 
equal-time structure factor S{ci2d) fl. In Fig. 1, we 
present representative scans of this type that show a 
peak in S'(q2D) the incipient antiferromagnetic or- 
dering wavevector. We fit these data to an isotropic 
2D Lorentzian, S{q2D) = S{Q)/{1 + q^j^^^), convolved 
with the instrumental resolution. As the temperature 
or the doping level is increased, the correlation length ^ 
decreases, which manifests itself as a broadening of the 
peaks. 

Figure 2 summarizes the behavior of ^(x,T) for sev- 
eral Ce concentrations. Over a wide temperature range, 
the spin correlations exhibit an exponential dependence 
on inverse temperature, ^{x,T) ~ exp [27rps(x)/T] with 
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FIG. 2: Semi-log plot of the magnetic correlation length (in 
units of the lattice constant a) versus inverse temperature (in 
units of J = 125 meV « 1450 K) for x = 0,0.10,0.14, and 
0.18. The data terminate just above the onset of Neel order at 
low temperature. Also displayed are previous data for a; = 
and X = 0.15 0,0 and Quantum Monte Carlo results for Eq. 
(1). 



spin stiffness Psix), indicative of an underlying ground 
state with 2D antiferromagnetic order We overlay 

the results of quantum Monte Carlo simulations for the 
randomly-diluted SLHAF, described by the Hamiltonian 



= J y^^piPjSj 

{i,3) 



^3 ' 



(1) 



where the sum is over NN sites, J is the antiferromag- 
netic Cu-O-Cu superexchange. Si is the spin-1/2 oper- 
ator at site i, Pi — 1 on magnetic sites, and = on 
nonmagnetic sites. The numerical method is the same 
as that employed previously Q. The data are plot- 
ted as ^/a, the ratio of the spin-spin correlation length 
to the NN Cu-Cu distance (a « 3.95 A), versus J/T, 
the ratio of the NN antiferromagnetic superexchange 
of Nd2Cu04 to the temperature. A previous estimate 
from comparison of £,(T) with theory for the SLHAF 
yielded J « 130 meV 0, and we find J ^ 125(4) mcV 
from comparison with our numerical results for a: = 0. 
Since the value of J has been established, the only ad- 
justable parameter in the comparison for the Ce-doped 
samples is the concentration of nonmagnetic sites in Eq. 
(1). We find good qualitative agreement between ex- 
periment and Monte Carlo if we fix this concentration 
to equal the nominal Ce concentration x. However, 
as shown in Fig. 2, the agreement becomes quantita- 



3 



tive if we allow the effective dilution to differ slightly: 
Xesix = 0.10) = 0.08(1), x^s{x = 0.14) = 0.12(1), and 
Xes{x = 0.18) = 0.21(2). Also shown are previous data 
for as-grown samples with a; = and x = 0.15 0,0 ■ The 
result for Nd2Cu04 agrees well with ours, and we esti- 
mate that Xcff{x — 0.15) — 0.16(2). The extent to which 
the quenched-disorder model Eq. (1) effectively describes 
the spin correlations of as-grown, non-superconducting 
NCCO is remarkable. The observed small differences be- 
tween nominal electron concentration and effective static 
dilution may be in part due to the uncertainty in the 
oxygen (and hence electron) concentration of our sam- 
ples. Another likely possibility is the presence of addi- 
tional quantum fluctuations which are not captured in 
the effective Hamiltonian Eq. (1), for example due to 
frustrating next-NN interactions. 

The spin stiffness of the disorder- free spin- 1/2 NN SL- 
HAF is known to be 2Trps/J ^ 1.131(3) An es- 

timate for the case of quenched disorder, Eq. (1), has 
been obtained in previous numerical work jl4||. Com- 
parison of the experimental results in Fig. 2 gives the 
effective values 27rp,/J = 0.71(2), 0.54(4), 0.40(7), and 
0.25(5), respectively, for x = 0.10, 0.14, 0.15, and 0.18. 
The value 0.40(7) for the x = 0.15 sample is somewhat 
smaller than the original estimate of 27rps/ J « 0.54 M 
which was based on a classical rather than a quantum |l4| 
picture for the doping dependence of the spin stiffness. 

We note that Ref. @ reports data for two as-grown 
x = 0.15 samples with onset of magnetic order at Tn — 
160 and 125 K. If we allow for a distribution of Neel tem- 
peratures to account for the rounding of the observed 
transition, we estimate mean values of Tn = 147(17) and 
115(9) K, respectively. This is consistent with how we 
have treated our own data. The second sample from Ref. 

also exhibits lower values of £,{T) (not shown in Fig. 
2), comparable to our x = 0.18 crystal. The removal of 
oxygens through the reduction process changes both the 
degree of disorder and the carrier concentration. This 
lowers the Neel temperature and, for x > 0.13, leads 
to a decrease of the antiferromagnetic volume fraction 
and to bulk superconductivity illj]. Growth at lower 
oxygen partial pressure mimics the reduction process, 
resulting in lower oxygen contents and higher effective 
doping levels. The primary focus of this Letter is on 
non-superconducting, as-grown NCCO, and therefore we 
were careful to grow our samples under high oxygen par- 
tial pressure, resulting in crystals with Tn at or close to 
the maximum value attainable. 

Figure 3(a) indicates that the Neel temperature Tjq{x) 
of our as-grown samples, as obtained from neutron 
diffraction, approximately follows a parabolic form, ex- 
trapolating to zero at a; « 0.21. As shown in Fig. 3(a), 
qualitatively similar behavior was found previously for 
reduced NCCO for which Tn{x) is lower The pri- 

mary effect of the reduction is an approximately rigid 
shift by |Aa;| « 0.03 due to the change in carrier density. 




FIG. 3: (a) Neel temperature doping dependence. Ce-doped 
samples exhibit some rounding in the transition due to oxy- 
gen and cerium inhomogeneities, which is represented by the 
error bars. Results from a previous study ^3| of reduced sam- 
ples are shown for comparison (triangles). The dashed curves 
are quadratic fits, and the dotted line indicates the behavior 
observed for LCZMO I3- The grey circles are the as-grown 
data shifted by |Aa::| = 0.03. (b) Doping dependence of the 
ordered copper moment per site, normalized by the value for 
as-grown Nd2Cu04. The grey circles are the as-grown data 
shifted by |Aa;| = 0.03. The squares are data obtained in a 
previous study and the dot-dashed line is a recent theo- 
retical Drediction [2(l I2I ^ . For comparison, Monte Carlo 
results for the randomly diluted NN spin-1/2 SLHAF (solid 
line) and the behavior observed for LCZMO (dotted line) 
[T^ are indicated as well. 



The concave doping dependence for NCCO is contrasted 
by the behavior exhibited by the randomly-diluted spin- 
1/2 SLHAF La2Cui_:r(Zn,Mg)^04 (LCZMO), for which 
Neel order extends up to the site percolation threshold, 
Xp « 0.41, and the 2D spin correlations are quantitatively 
described by Eq. (1) up to Xp [l^ . 

In Fig. 3(b), we plot the ordered moment for sev- 
eral samples, as obtained from order parameter mea- 
surements at low temperature. The values are normal- 
ized to that of Nd2Cu04. Consistent with the behav- 
ior of Tf^{x), the ordered moment of as-grown NCCO 
approaches zero much more rapidly than in the case of 
random dilution: at the rate 1 — Mst(a;)/Mst(0) w 3.5x 
which is approximately twice that of the diluted spin-1/2 
SLHAF for x < 0.20. On the other hand, the magnetic 
phase extends much further than for LSCO and Li- 
doped La2Cu04 for which Neel order is destroyed 
above x w 0.02 and x « 0.03, respectively. The ordered 
moment of reduced samples is lowered further. We note 
that these data are again normalized by the value for as- 
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grown Nd2Cu04, and that a rigid shift of |Aa;| = 0.03 
leads to a good agreement of the two data sets. Our 
result is consistent with previous experimental work at 
lower concentrations [Toj . 

If viewed instantaneously, so as to mitigate the ef- 
fects due to the increased itinerancy upon doping, the 
magnetism of NCCO should indeed resemble a system 
with quenched random site dilution. What is surpris- 
ing is that for ^{x, T) the correspondence with Eq. (1) is 
nearly quantitative. Since the charge carriers in as-grown 
NCCO are itinerant, even well below Tn [l^l, deviations 
from the static model Eq. (1) may be expected to emerge 
as we view the two systems on different time scales. In- 
deed, the ordered moment, a measure of the strength 
of the magnetic order at infinitely long times, decreases 
much more rapidly in the case of NCCO than for the SL- 
HAF. We note that slight deviations, possibly due to the 
presence of a frustrating next-NN exchange, are already 
observable for the ordered moment of LCZMO [Fig. 
3(b)], and that LCZMO appears to be in close proximity 
to a new quantum critical point in an extended param- 
eter space 10. In the case of LCZMO, these pertur- 
bations are not quite strong enough to shift the critical 
point away from the geometric site-percolation threshold. 
Although the experimentally accessible doping range is 
limited, our results for NCCO are consistent with the ex- 
istence of a quantum critical point below the geometric 
site-percolation threshold. 

We note that recent calculations for the one-band Hub- 
bard model are in semi-quantitative agreement with our 
data for the sran correlations [20I Esf and ordered mo- 
ment 113,111123. The latter result is shown in Fig. 3(b). 
These theories were motivated by photoemission mea- 
surements of the evolution of the Fermi surface of reduced 
Nd2-xCe:^Cu04±5 (0 < X < 0.15) H^. Their proponents 
assert that t-J models are intrinsically limited in describ- 
ing the electron-doped cuprates and that a proper de- 
scription requires t-U models with a decreasing value of 
U with increasing electron concentration [2iJ, |2l|, |22f . 

In summary, we have presented a comprehensive neu- 
tron scattering study of the magnetic properties of non- 
superconducting Nd2_a;Cea;Cu04±5. Although the evo- 
lution of the instantaneous spin-correlations may effec- 
tively be described by a site-dilution model, a proper 
description of the magnetic degrees of freedom must nec- 
essarily include the effects of electron itinerancy. While 
the achievements of recent theoretical treatments of the 
Hubbard model seem very promising, it is not clear at 
present if this model exhibits a superconducting ground 
state. Phonon anomalies have been observed in both 
electron and hole doped high-temperature superconduc- 
tors ll^, and a complete description of these materials 
may require the additional consideration of the electron- 
lattice coupling. Regardless, our data should serve as 
a benchmark for tests of still-emerging theories for the 
high-Tc phase diagram. 
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